Flow at the SPS and RHIC as a Quark-Gluon Plasma Signature 
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Radial and elliptic flow in non-central heavy ion collisions can constrain the effective Equation of 
State(EoS) of the excited nuclear matter. To this end, a model combining relativistic hydrodynamics 
and a hadronic transport code(RQMD Jl|) is developed. For an EoS with a first order phase 
transition, the model reproduces both the radial and elliptic flow data at the SPS. With the EoS fixed 
from SPS data, we quantify predictions at RHIC where the Quark Gluon Plasma(QGP) pressure is 
expected to drive additional radial and elliptic flow. Currently, the strong elliptic flow observed in 
the first RHIC measurements does not conclusively signal this nascent QGP pressure. Additional 
measurements are suggested to pin down the EoS. 

1. By colliding heavy nuclei at the SPS and RHIC 
accelerating facilities, physicists hope to excite hadronic 
matter into a new phase consisting of deconfined quarks 
and gluons - the Quark Gluon Plasma(QGP) |1. After 
the collision, the produced particles move collectively or 
flow and this flow may quantify the effective Equation 
of State(EoS) of the matter. In central PbPb collisions 
at the SPS, a strong radial flow is observed ||. The mat- 
ter develops a collective transverse velocity approaching 
(l/2)c. In non-central collisions, a radial and an elliptic 
flow are observed Since in non-central collisions 

the initial nucleus-nucleus overlap region has an elliptic 
shape, the initial pressure gradient is larger along the 
impact parameter and the matter moves preferentially in 
this direction jjj. 

The phase transition to the QGP influences both the 
radial and elliptic flows. QCD lattice simulations show 
an approximately 1st order phase transition Over a 
wide range of energy densities e = .5 — 1.4 GeV/ fm 3 , the 
temperature and pressure are nearly constant. Over this 
range then, the ratio of pressure to energy density p/e, 
decreases and reaches a minimum at a particular energy 
density known as the softest point, e sp « 1.4 GeV/ fm 
H . When the initial energy density is close to e sp , the 
small pressure (relative to e) cannot effectively acceler- 
ate the matter. However, when the initial energy den- 
sity is well above e sp , p/e approaches 1/3, and the larger 
pressure drives collective motion ||,[l0|]. At a time of 
~ 1/m/c, the energy densities at the SPS(y / i 7VAr = 
17 GeV) and RHIC (\^s NN = 130 GeV) are very approx- 
imately 4 and 7 GeV/ fm 3 respectively fi"l"l , |l2"f . Based on 
these experimental estimates, the hard QGP phase is ex- 
pected to live significantly longer at RHIC than at the 
SPS. The final flows of the produced particles should re- 
flect this difference. In this paper we pose the question: 
Can both the radial and elliptic flow at the SPS and 
RHIC be described by a single effective EoS? 

Since the various hadron species have different elastic 
cross sections, they freezeout (or decouple) from the hot 
fireball at different times [jl^|. Because flow builds up 
over time, it is essential to model this differential freeze- 
out. It was ignored in previous hydrodynamic simula- 
tions of non-central heavy ion collisions and elliptic flow 



was over-predicted flow by a factor of two [^4||Hj]. 

2. The Hydro to Hadrons(H2H) model will be de- 
scribed in detail elsewhere |l(| . Other authors have pre- 
viously constructed a similar model for central collisions 
p"7| . The model evolves the QGP and mixed phases 
as a relativistic fluid, but switches to a hadronic cas- 
cade (RQMDv2.4 0) at the be ginning of the hadronic 
phase to model differential freezeout. The computer 
code consists of three distinct components. Assuming 
Bjorken scaling, the first component solves the equations 
of relativistic hydrodynamics in the transverse plane Q| 
and constructs a switching surface at a temperature, 
Tswitch = 160 MeV. The second component generates 
hadron on the switching surface using the Cooper-Frye 
formula Jl8| with a theta function rejecting backward 
going particles [^9 20 1. Finally, the third component 
(RQMD) sequentially re-scatters the generated hadrons 
until freezeout. 




e (GeV/fm ) 

FIG. 1. The pressure versus the energy density(e) for dif- 
ferent EoSs (see text). EoSs with Latent Heats 0.4 GeV/ fm 3 , 
0.8GeV/fm 3 ,... are labeled as LH4, LH8,...etc. 

For the hydrodynamic evolution, a family of EoSs 
was constructed with an adjustable Latent Heat(LH) 
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FIG. 2. The transverse mass slope(T s ; ope ) as a function of 
the total charged particle multiplicity in PbPb collisions at an 
impact parameter of b=6 fm (see also [ fu)||l7| ) . For consistency 
with the elliptic study in Fig. |§| we show b=6 fm although the 
NA49 data points fell are for the 5% most central events, or 
b<3.5 fm. For all EoSs at the SPS, the proton slope param- 
eters at b=6fm are ~ 7 MeV smaller than at b=0fm, as for 
the b=0 LH8 curve. The difference is negligible for n~~ . 

(see Fig. [j]). LHoo is considered as a limiting case, 
mimicking non-equilibrium phenomena [^2| . The hadron 
phase exists up to a critical temperature of T c — 
165 MeV, and consists of ideal gas mixture of the me- 
son pseudoscalar and vector nonets and the baryon octet 
and decuplet. The hadron phase is followed by a mixed 
phase with a specified LH, which is finally followed by 
a QGP phase with C 2 = 1/3. In addition, a Resonance 
Gas(RG) EoS was constructed with a constant speed of 
sound above the hadron phase. 

3. Radial flow is quantified experimentally by slope 
parameters, T s i ope ; the momentum spectrum of each par- 
ticle is fit to the form dN/dM T dy\ y=0 = Ce' UT/T "°' e 
where M T = P T + m 2 . T s i ope incorporates random ther- 
mal motion and the collective transverse velocity. 

In Fig. ||, the pion and proton slope parameters are 
plotted as functions of the total charged particle multi- 
plicity in the collision. Look first at the leftmost points 
at SPS multiplicities and compare the model and exper- 
imental slopes: The proton slope data favor a relatively 
hard EoS - LH8 or harder. A direct comparison of the 
model to published spectra ^3j supports this claim [jl6| . 
A RG EoS can also reproduce the proton flow. A simi- 
lar analysis of elliptic flow (shown and quantified below) 
favors a relatively soft EoS - LH8 or softer. With some 
caveats, LH8 represents a happy middle which can repro- 
duce both the radial and elliptic flow at the SPS. 

Look now at the energy/multiplicity dependence of the 
slopes. For all EoSs, T s i ope increases with the collision 
energy flfl,0]. For a soft EoS (e.g. LHoo) the increase is 



small and for a hard EoS (e.g. LH8) the increase is large. 
At RHIC multiplicities, the difference between the slope 
parameters is large and easily experimentally observable. 

4. Elliptic flow is quantified experimentally by the el- 
liptic flow parameter, V2 = (cos(2$)); here $ is the an- 
gle around the beam measured relative to the impact 
parameter and ( ) denotes an average over the single 
particle distribution, d p^ d<i , ■ V2{Pt) is found by hold- 
ing Pt constant in while averaging cos(2<I>) over d p^ di - 
V2 measures the response of the fireball to the the spa- 
tial deformation of the overlap region, which is usually 
quantified in a Glauber model |2J] by the eccentricity 
e = ((y 2 — x 2 )) / ((x 2 + y 2 )) . Since the response^) is pro- 
portional to the driving force(e), the ratio v-ij t is used to 
compare different impact parameters and nuclei [ p5|j26| . 

In Fig. |^(a), the number elliptic flow (t> 2 ) is plotted as 
a function of charged particle multiplicity at an impact 
parameter of 6 fm. Before studying the energy depen- 
dence, look at the magnitude of the elliptic flow at the 
SPS. For LH8, the stars show the pion V2 when the mat- 
ter is evolved as a fluid until a decoupling temperature 
of Tf = 120 MeV; they illustrate the excessive elliptic 
flow typical of pure hydrodynamics. Once a cascade is 
included, LH8 (the squares) is only ps 20% above the 
data - a substantial improvement. Typically in hydro- 
dynamic calculations, the freezeout temperature Tf is ad- 
justed to fit the proton Pt spectrum. However, protons 
are driven by a pion "wind" and decouple from the fire- 
ball 5 fm/c after the pions on average. This pion wind 
accounts for the strong proton flow at the SPS and is 
not described by ideal hydrodynamics [|l7|Jl6| . In order 
to match the observed proton flow, hydrodynamic cal- 
culations must decouple at low freezeout temperatures, 
Tf rz w 120 MeV/c. This low temperature has two con- 
sequences for elliptic flow: First the reduction of elliptic 
flow due to resonance decays is small ~ 15%, compared 
to ss 30% in the H2H model. Second, compared to a 
cascade, the hydrodynamics generates twice as much el- 
liptic flow during the late cool hadronic stages of the 
evolution. By including the pion "wind" , and more gen- 
erally by decoupling differentially, we can simultaneously 
describe the radial and elliptic flow data at the SPS. 

The energy dependence of «2 is the central issue. As 
seen in Fig. |3|, the H2H model predicts an increase in 
elliptic flow by a factor w 1.4 and is in reasonable agree- 
ment with SPS and RHIC flow data. This result was 
presented prior to the publication of RHIC data J2(| . In 
contrast, UrQMD, a hadronic cascade based on string 
dynamics, predicts a decrease by a factor of w 2 p7) . 
This is because the UrQMD string model has a super- 
soft EoS at high energies p8[ . For pure hydrodynamics 
as illustrated by the stars, V2 is approximately constant 
|4j (but see (15)). For HIJING |§, a model which con- 
siders only the initial parton collisions, V2 is ~ [ Btj . 
The first RHIC data clearly contradict these models. 

The increase in V2 is now used to constrain the EoS of 
the excited matter. The QCD phase diagram has two 
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FIG. 3. (a) The number elliptic flow parameter V2 as a 
function of the charged particle multiplicity in PbPb colli- 
sions at a an impact parameter of b=6fm. At the SPS, the 
NA49 V2 data point is extrapolated to b=6fm using Fig. 3 
in ||. At RHIC, the STAR V2 data point is extrapolated to 
N ch /N^ ax = 0.545 (b=6fm in AuAu) using Fig. 3 in [|. The 
comparison to data is a little unfair: For the model, V2 is cal- 
culated using all pions in PbPb collisions. For the NA49 data, 
V2 is measured using only tt~ in PbPb (a -3% correction to 
the model). For the STAR data, 1)2 is measured using charged 
hadrons in AuAu (a +5% correction to the model). 

distinguishing features. It is soft at low energy densities 
and subsequently hard at high energies. A RG EoS (the 
open squares) has no softness and the elliptic flow is 
clearly too strong both at the SPS and RHIC. The entire 
family of EoSs, LH8 through LHoo, reproduces the ellip- 
tic flow data in both energy regimes. Counter-intuitively, 
as the latent heat is increased, v 2 first decreases and then 
increases. In the final count, LH8 and LHoo have roughly 
the same v 2 . However, they develop the v 2 in differ- 
ent ways. For LH8, the EoS shifts from hard to soft 
and the early pressure starts an early elliptic expansion. 
For LHoo, the EoS is just soft and the elliptic expan- 
sion stalls. However because the expansion is stalled, the 
LHoo collision lifetime^ 13 fm/c at RHIC) is signifi- 
cantly longer than the LH8 lifetime(« 9 fm/c at RHIC) 
P). Over the long LHoo lifetime, v 2 H °° steadily grows 
and is finally comparable to v 2 Ha - As the latent heat is 
increased from LH8 to LH16, the EoS becomes becomes 
softer and v 2 at first decreases. However, as the EoS is 
made softer still, the lifetime increases and v 2 rises again. 

5. Impact Parameter Dependence. In Fig. ^, v 2 for LH8 
as a function of the number participants (N p ) is compared 
to data. Different EoSs show a similar participant (or b) 
dependence. The agreement is good at RHIC where the 
multiplicity is high. For ideal hydrodynamics, v 2 oc e oc 
(N™ ax - N p ) 0. In the low density limit, since the 
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FIG. 4. V2 versus the number of participants(iVp) relative 
to the maximum . The model and the NA49 V2 values Q 
at the SPS are for n~ . The NA49 data are mapped from 
b to participants using |3lJ. The model and the STAR V2 
values H at RHIC are for charged particles. The model does 
not include weak decays. The number of charged particles is 
assumed proportional to N p . 



response is proportional to the number of collisions, 
v 2 oc oc {N p nax ~ N p )N p . Therefore, v 2 has a differ- 
ent N p (oT b) dependence in the hydrodynamic and low 
density limits |^,|6|. At RHIC, except in very periph- 
eral collisions, the N p dependence is clearly linear and 
strongly supports the hydrodynamic limit |^|. At the 
SPS, the N p dependence may not be clearly linear, but 
it also does not follow the low density limit. Two-pion 
correlations, may change the data analysis jj2|, reduce v 2 
in the periphery and improve the low density agreement. 

Finally in Fig. ||, v 2 is studied both as a function of 
transverse momentum and impact parameter. For both 
LH8 and LHoo, the calculation produces too much ellip- 
tic flow in peripheral collisions(45-85%), and too little 
elliptic flow in the most central collisions(0-ll%). The 
Pt dependence of v 2 also clarifies the difference between 
and LH8 and LHoo: LHoo, a super soft EoS, generates 
elliptic flow only at low momentum while LH8, a hard 
EoS, generates elliptic at high momentum. 

6. Summary and Discussion. By incorporating differ- 
ential freezeout, the Hydro to Hadrons (H2H) model si- 
multaneously reproduces the radial and elliptic flow at 
the SPS and RHIC. At the SPS, the radial flow de- 
mands an EoS with a latent heat LH > 0.8 GeV/ fm 3 , 
while elliptic flow demands an EoS with a latent heat 
LH < 0.8 GeV/ fm 3 . Further, in contrast to string and 
collision-less parton models, the increase in v 2 is natu- 
rally explained using hydrodynamics. This challenges the 
prevailing view |p§ that the SPS is in the low density 
regime and that the increase in v 2 represents a transition 
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FIG. 5. Elliptic flow of charged pions as a function of Pt 
and centrality for AuAu collisions at RHIC. For centrality se- 
lections The percentages shown 0-11%, 11-45% and 45-85%, 
indicate the fraction of the total geometric cross section for 
three centrality selections fm<b<4.2 fm, 4.2 fm<b<8.4 fm 
and 8.4 fm<b<11.6 fm. The preliminary data points were 
presented in fi3| , fj4"[ . The model curves were found by pa- 
rameterizing the model data points and averaging over the 
specified impact range with the geometric weight, 2-ubdb. 

to the hydrodynamic regime. However, the increase in V2 
does not uniquely signal the asymptotic QGP pressure. 
Indeed, at RHIC collision energies, a very soft EoS can 
have the same v-i as an EoS with a well developed QGP 
phase. This EoS is not academic since softness can mimic 
non-equilibrium phenomena p2[ . To reveal the underly- 
ing EoS and the burgeoning QGP pressure, the collision 
energy should be scanned from the SPS to RHIC. If the 
prevailing low density view of the SPS is correct, a tran- 
sition in the b dependence of elliptic flow should be ob- 
served over the energy range p5| , |26| . In addition, for 
different EoSs, depends differently on collision energy 
and transverse momentum (Fig. || and Fig. |^). Taken 
with the radial flow (Fig. this experimental informa- 
tion would help settle the EoS of hot hadronic matter. 
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Note Added. After the submission of this work, the 
STAR and PHENIX collaborations reported proton at 
anti-proton spectra p3[ |. The preliminary spectra favor 
LH8-LH16 and disfavor LHoo fell. 
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